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Simulation  of  Repetitively-Pulsed  Laser  Irradiation 
of  Graphite-Epoxy  Composite 


1.  INTRODUCTION 

For  their  investigations  of  the  rapid  heating  of  graphite  epoxy 
composites,  Griffis,  Masumura  and  Chang  (GMC)  [1]  developed  a  computer  code 
to  calculate  the  response  of  a  composite  to  CW  laser  irradiations  with 
fluences  of  the  order  of  kilowatts  per  square  centimeter.  Data  are  also 
available  for  the  repetitively-pulsed  laser  irradiations  of  the  same 
composite  material.  These  Irradiations  involve  much  higher  peak  powers,  of 
the  order  of  megawatts  per  square  centimeter,  typically  with  100  oulses  of 
10  us  duration  at  the  rate  of  100  pulses  per  second.  The  average  energy 
deposited  (ICO  X  lC'1’  s  X  106  W/cm2  «  103  J/cn2)  Is  about  the  same 
as  In  the  GMC  work. 

In  this  report  I  describe  the  modifications  that  I  made  to  the  GMC  code 
so  as  to  handle  the  larger  pulsed  fluences.  When  attempting  to  use  the  GMC 
code  at  higher  fluerxes,  reveral  problems  arose.  The  code  had  to  be 
modified  to  handle  pulsed  Irradiations,  but  this  did  not.  In  itself,  disturb 
the  numerical  procedures  used  In  the  code.  The  next  section  describes  the 
new  time-mesh  structure,  which  allows  a  repeated  heat-on/heat-off 
Irradiation.  The  three  order-of-maghltude  Increase  In  the  peak  fluence  did 
generate  numerical  problems.  Oue  to  the  low  thermal  conductivity  of  the 
composite  material,  enormous  thermal  gradients  are  created  at  the  surface  of 
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the  material— gradients  of  the  order  of  10  -1G  C/cm.  This  problem  was 
circumvented  by  replacing  the  grid  sturcture  for  the  depth  variable  with 
another  structure  that  has  a  very  fine  mesh  size  at  the  surface  and  that 
gradually  coarsens  as  one  goes  deeper  Into  the  material.  This  mesh  Is 
described  In  Section  ?. 
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In  the  GMC  model  of  the  composite  response  at  high  temperature,  the 
material  (graphite  and  residual  epoxy)  decomposes  at  3316 °C .  Once  the  front 
surface  reaches  this  temperature,  the  surface  begins  to  recede.  CMC 
modified  the  equations  describing  the  surface  conditions. in  order  to 
simulate  this  recession.  At  the  higher  fluences  used  with  the  repetitively- 
pulsed  irradiations,  their  procedure  becomes  numerically  unstable.'  I 
describe  in  Section  4  a  modified  method  for  handling  the  recession. 

Section  5  describes  the  changes  made  in  the  handling  of  the  thermal 
properties,  which  has  the  effect  of  both  improving  the  numerical  convergence 
of  the  calculations  and  improving  their  accuracy.  The  following  section 
reports  a  test  of  the  code  that  compares  the  program’s  calculation  of  the 
time  for  a  complete  burnthrough  with  the  time  derived  from  simple  heat 
capacity  estimates. 

in  order  to  verify  that  none  of  the  changes  outlined  above  made  any 
substantial  change  In  the  low  fluence,  CW  results.  Section  7  compares  one  of 
the  calculations  made  by  GMC  to  the  same  calculation  made  with  the  modified 
code.  The  next  section  contains  a  sample  calulation  for  a  pulsed 
irradiation.  Section  9  then  presents  a  detailed  comparison  of  this  code's 
results  with  an  experimental  pulsed  irradiation. 

Throughout  this  report  we  follow  the  notation  of  Griffis,  Masumura  and 
Chang.  Rather  than  repeat  equations  that  appear  in  their  report,  we  will 
refer  to  them  directly  by  using  their  equation  numbers  prefaced  by  a  "G". 

We  also  follow  GMC  In  using  the  expressions  for  surface  losses  due  to 
convection  and  reradiation  that  are  given  by  Hobbs,  et  al  [2].  The' 
concluding  section  will  emphasize  the  need  for  better  measurements  of  the 
thermal  properties  of  the  composite. 

2.  THE  TIME  EVOLUTION  MESH 

Ue  describe  the  structure  of  the  a  pulsed  irradiation  by  three  time 
periods:  tjf  the  length  of  an  Individual  User  pulse;  t»,  the  time 
between  pulses;  and  t^,  the  length  of  time  the  coce  is  allowed  to  run  at 
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the  end  of  the  irradiation.  Thus  the  total  time  simulated  by  the  code  is 

t  .  np  (tx  +  t2)  ♦  t3  (1 

where  np  is  the  number  of  pulses.  Typically,  an  experimental  run  will 
consist  of  100  pulses  at  100  pps,  with  1C  ps  pulse  lengths,  for  which  case 
we  have  tj*10  ps,  t2=9990  ws,  and  np=10C.  Having  set  these  intervals, 
we  break  them  into  n^,  n^,  and  n^  subintervals, 


*s  *  *sl  *  ts2  *  •**  *  tsn$* 

so  that  there  are  a  total  of  n^  time  intervals,  where 


nt  *  np  (nl  *  n2}  +  n3  . 


The  code  is  most  unstable,  numerically,  at  the  beginning  of  a  pulse--the 
thermal  gradients  are  largest  then  and  the  difficulties  associated  with  the 
abrupt  changes  in  the  composite's  thermal  properties  [1]  are  worst. 
Consequently,  the  code  requires  small  time  intervals  at  the  beginning  of  a 
pulse.  Later,  after  the  heat  has  flowed  into  the  interior  of  the  composite 
and  after  the  temperature  has  risen  above  the  region  (2C0-5CC*C)  where  the 
thermal  properties  are  changing  rapidly,  larger  time  intervals  can  be  used. 
We  satisfiy  the  requirements  by  setting  the  time  steps  for  the  pulse  period 
as 


expii  a'j)  -  exp({1-l)  aj) 


with  similar  expressions  for  the  otheh  time  Intervals.  The  values  of  a^, 
a2,  and  a^  can  be  set  as  desired;  typical  values  are  a^  «  0.07, 
a2  «  C.10,  and  a ^  >  C.10,  in  which  case.  If  t^*10  ps,  n^.lCC, 
t2«5.99  ms,  and  n2»10G,  the  first  tint  step  is  0.66  ns  long,  the  ICOth 
is  0.68  us,  the  101st  Is  48  ns,  and  the  200th  Is  0.95  ms.  The  form  (4)  is 
flexible,  allowing  great  disparity, in  the  time  steps,  as  above,  and  allowing 
equal  steps  as  the  a^  go  to  zero. 
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In  the  course  of  its  calculations,  if  the  code  runs  into  numerical 
difficulty  during  a  particular  time  step,  such  as  the  non-convergence  of  one 
of  the  iteration  procedures,  then  the  size  of  that  time  step  is  halved  and. 
the  procedure  restarted.  Once  the  whole  time  step  has  been  worked  through, 
the  program  reverts  to  the  original  time  step  sequence.  Thi*  self  adjusting 
feature  makes,  the  code  more  flexible,  so  that  larger  and  more  efficient  time 
steps  can  be  used  most  of  the  time;  shorter  time  steps  are  used  only  when 
necessary. 


3.  THE  MESH  FGR  PENETRATION  DEPTH 

When  a  several  megawatt  irradiation  first  impinges  on  the  composite 
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material,  thermal  gradients  of  the  order  of  10  -1C  C/cm  are  seen.  In 
order  that  the  calculated  temperatures  not  vary  too  much  from  gird  point  to 
adjacent  gird  ooint,  it  is  necessary  that  the  mesh  be  quite  fine  near  the 
surface.  To  accomplish  this  we  again  use  the  exponential  form 

Xq  *  C 

,x.j  *  L  [exp(i  b)  -  1)  /  (expln^  b)  -  1]  ,  1=1,2. ,.nm  (5 

where  l  is  the  thickness  of  the  material  and  n  is  the  initial  number  of 

m 

mesh  points.  Typically,  we  use  a  value  of  b*C.035,  so  that  for  a  thickness 

Q 

of  0.254  cm  the  first  grid  slab  has  a  thickness  of  4X10"  cm  and  the  final 

_*3 

slab  a  thickness  of  5X^0  ~  cm. 

4.  NEW  ITERATI'T  METHOO  DURING  ABLATION 

The  equations  governing  the  response  of  the  composite  change  once  the 
front  surface  of  the  composite  reaches  the  ablation  temperature.  The 
equation  [1]  relating  the  surface  temperatures  and  the  recession  velocity  is 
Eq.  (G15),  which  we  reproduce  here 


Tlj  +  Tljn  ■  T2j  ■  T2j+1  = 

=  cf  [Zzz  k2  -  P  cp  k  Vm  z 2  -  Qf  (dk/dT)  z2}  k"~  (6 

Qf  =  C-'pHsVm  ,  (7 

where  Q  is  the  absorbed  power  flux  less  the  surface  losses.  (The  second 
Cy  factor  in  Eq.  (6)  was  inadvertently  left  out  of  Eq.  (15)  of  Ref.  ( 1 ) . ) 

Griffis,  Masumura  and  Chang  set  the  frorit  surface  temperature  during 
ablation  to  be 


They  then  solve  the  coupled  equations  for  the  new  set  of  terrperatures 
[T.Jn,  i=2,...,nm],  and  then  use  Eq.  CGI 5)  (with  T^.  =  T^+1  = 

T$)  to  determine  V^.  They  then  iterate  using  the  new  Vm  and  T^+j 
values.  This  process  is  repeated  until  successive  values  of  Vm  are  within 
2  percent  of  each  other. 


When  I  attempted  to  use  their  iteration  scheme  with  higher  fluences 

2  1 
(2  MW/cm  ),  the  code  became  unstable,  so  I  created  a  modified  interation 

method  based  on  the  same  equations.  Instead  of  Eq.  (8 ) ,  I  used  Eq.  (7)  as 

one  of  the  nm  equations  ,1  inking  the  temperatures  at  the  nm  grid  points. 

After  the  new  set 

the  surface 

prescribed 

time  step  is 


The  procedure  is  started  by  guessing  at  a  value  for  \ 
of  temperatures  [T^.+j,  1*1,2, . . . ,nm3  is  calculated, 
temperature  T^^  is  checked  to  see  if  it  is  within  a| 
neighborhood  of  T$.,  If  so,  the  calculation  for  that 
finished.  If  not,  the  new  temperatures  are  scaled 


T’lj*l  *  (Ts  1  Tlj+1}  Tij*l 


and  a  new  guess  Is  made  as  to  the  value  of  Vm<  This 
until  a  given  Vffl  yields  T^+j  within  the  prescribed 
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procedure  is  repeated 


neighborhood  of  T  , 


> 


i 


vTI 


5.  NEW  f-iETHCD  FOR  REPRESENTING  THERMAL  PROPERTIES 

The  thermal  properties  we  use  in  these  calculations  are  those  given  in 
Menousek  and  Monin  [3].  The  density,  specific  heat,  and  thermal 
conductivity  are  represented  by  a  series  of  ramp  and  step  functions;  the  low 
energy  portions  of  these  properties  are  pictured  in  the  following  figures. 

Figure  1  depicts  the  Menousek  and  Monin  characterization  of  the  specific 
heat  of  the  graphite  epoxy  as  it  is  heated  from  G*C  to  900*0.  The  large 
changes  in  the  specific  heat  between  340*0  and  510*0  reflect  the  chemical 
changes  that  the  epoxy  undergoes  between  these  temperatures.  What  remains 
above  510*0  are  the  graphite  fibers  ar<d,  presumably,  some  residue  of  the 
pyrolyzed  epoxy.  Although  it  is  not  completely  shown  in  Figure  1,  there  is 
.a  gentle  rise  of  the  specific  heat  between  510*0  and  3316*0  where  the 
graphite  fibers  and  epoxy  residue  sublimate.  We  have  used  the  Menousek  and 
Monin  values  for  the  specific  heat,  with  one  modification.  Once  the 
composite  has  been  heated  above  510*0  the  epoxy  is  permanently  lost,  so 
thereafter  the  specific  heat  is  represented  by  the  line  from  510*0  to  3316*0 
extended  down  to  0*0;  the  extension  is  shown  by  the  large-dash  curve  in 
Fig.  (1).  Also  shown  in  this  figure  is  the  characterization  used  by  GMC, 
which  consists  of  the  Menousek  and  Monin  curves  with  the  steeper  ramps 
replaced  by  gentler  ones.  We  will  discuss  this  point  further  after 
displaying  the  other  properties. 

Figure  2  shows  the  Menousek  and  Monin  version  of  the  density  up  to 
90C*C;  It  is  constant  above  510*0.  We  again  use  their  version  upon  first 
heating  the  composite.  Once  the  material  has  been/heated  above  51C*C  we 
assume  a  constant  value  of  the  density  of  1.084  g/cm  .  The  slightly 
modified  GMC  version  for  the  density  is  also  shown  in  the  figure.  The  same 
discussion  applies  to  Figure  3,  showing  the  thermal  conductivity,  with  one 
exception.  The  Menousek  and  Monin  curve  Is  a  straight  line  dropping  from  a 
value  of  1.452  W/cm*C  at  10*0  to  a  value  0.173  W/cro*C  at  £38*0,  after  which 
It  is  constant  at  that  value.  For  convenience,  we  shifted  the  position  at 
which  the  curve  turns  from  538*0  to  51G*C,  the  same  temperature  at  which  the 
ether  properties  break.  Given  the  uncertainties  with  which  these  properties 
are  known  [3],  this  seems  a  negligible  change.  All  of  the  properties  that 


Fig.  1  —  Specific  Heat  of  the  graphite  epoxy.  Shown  1$  the  original 
Menousek  and  Monin  version,  which  Is  used  in  the  current  study,  and  the 
slightly  modified  version  used  by  GMC.  Also  shown  is  the  form  used  in  thl 
study  when  the  sample  has  been  heated  enough  to  pyrolyt*  the  epoxy. 


DENSITY  [g/cm3] 


have  been  presented  graphically  in  this  section  are  given  in  tabular  form  in 
Appendix  A. 


0  200  400  600  800 
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Fig.  2  --  The  density  of  the  graphite  epoxy..  Shown  is  the  original  Menousek 
and  Monin  version,  which  is  used  in  the  current  study,  and  the  slightly 
modified  version  used  by  CMC.  Also  shown  Is  the  form  used  in  this  study 
when  the  sample  has  been  heated  enough  to  pyrolyre  the  epoxy. 


THERMAL  CONDUCTIVITY  [W/cm°CJ 


Fig.  3  —  The  thermal  conductivity  of  the  graphite  epoxy.  Shown  Is  the 
version  used  In  this  study  and  the  version  used  by  GMC,  both  of  which  are 
slightly  modified  forms  of  the  original  Menousek  and  Honin  version.  Also 
shown  Is  the  form  used  In  this  study  when  thfe  sample  has  been  heated  enough 
to  pyrolyze  the  epoxy. 


We  conclude  this  section  on  thermal  properties  by  describing  an  improved 

method  of  utilizing  these  properties  within  the; code.  The  original  GMC  code 

supplies  tnermal  properties  in  the  following  manner  auring  a  given  time 

interval:  For  each  grid  step  of  the  spatial  coordinate,  there  is  an  initial 

temperature  T.  and  an  estimated  final  temperature  T_.  The  GMC  code 
J  .  1 

averages  these  two  temperatures  and  then  supplies  the  thermal  properties 

corresponding  to  the  average  temperature  by  interpolating  in  Table  A3.  A 

more  accurate  method  of  supplying  the  thermal  properties  is  to  provide  the 

average  thermal  properties  over  the  temperature  range  1\  to  T^.  As  an 

example,  we  will  use 


rather  than  using  the  thermal  conductivity  at  t<-?e  average  temperature 


(yv/2- 

Because  T 2  varies  with  each  Iteration,  If  one  of  the  temperatures  lies, 
near  one  of  the  snarper  boundaries  shown  in  Figures  1-3,  then  minor  shifts 
in  the  temperature  can  make  significant  shifts  n  the  reported  thermal 
properties.  This  effect  tends  to  destablize  th*.  convergence  of  the 
iteration  procedure. 

In  Figure  4  we  snow  an  example  of  the  benefit  of*  the  improved  method; 
the  specific  feat  is  plotted  as  a  function  of  T2  while  T.  is  fixed  at 
4SG*C.  In  this  example,  there  is  as  much  as  a  ten  percent  error  in  Cp  at 
T2»530*C  ( (T2^Tj )/2  •  510*C).  (Even  larger  errors  occur  in  this 
example  when  T,.=56G*Cj  but  usually  the  time  steps  are  picked  so  that  this 
large  a  change  in  temperature  would  not  arise.)  In  the  same  vein,  the 
derivative  of  the  diffuslvity,  which  appears  in  fcq.  (6),  has  discontinuities 
as  one  crosses  some  of  the  boundary  points.  The  GMC  method  of  supplying  the 
average  value  of  this  quantity  also  yields  discontinuities.  These 
discontinuities  can  lead  to  minor  instablitles  In  the  Iteration  process, 
requiring  a  greater  number  of  Iterations  and  producing  less  accuracy  in  the 
final  results.  The  new  method  described  in  this  section  provides  a  smoother 
and  a  more  realistic  averaging  process. 
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AVERAGE  SPECIFIC  HEAT  (490°C  -  T)  [J/°C  g] 


f  i  i 

- CURRENT  STUDY 

- GRIFFIS,  MASUMURA  &  CHANG 


Fig.  4  --An  example  of  the  specific  heat  as  supplied  by  the  current  study 
and  by  &MC.  The  text  above  describes  the  difference  between  the  two  methods. 


6.  SIMPLE  HEAT  CAPACITY  TEST 


If  the  specific  heat  anc  density  of  this  mcoel  material  were  constant 
and  7,  respectively ),  then  the  tine  that  is  needed  to  ccmoletely  burn 
through  a  sample  of  thickness  d  wculc  be 

a  o  (  1T~  iT  *  H  ) 

.  P  s  ,,, 


where  C  is  the  power  of  the  incoming  radiation  in  W  cm  ,  a  is  the 
absorption  coefficient,  represents  the  surface  losses,  aT  is  the 
temperature  rise  from  ambient  to  sublimation  of  the  graphite,  and  Hs  is 
the  heat  of  sublimation  of  the  graphite  fibers  and  epoxy  residue. 

The  average  velocity  of  recession  is 

Vave  ■  d  '  *8  .  112 

and  the  maximum  velocity  of  recession  is 

vr,aax  *  (  a  C  *  CL)  /  o  \  •  (13 

which  is  the  value  reached  just  before  bumthrough. 

In  Figure  5  we  show  the  surface  recession  velocity  as  <i  function  of  the 
time  of  irradiation  for  the  following  conditions: 

d  .  C.254  cm 

C  «  2.2  kU 
a-0.92 


0  hs  .  1.084  g/cm-  X  43  kJ/g 


■  46.6  kJ/cm 


«_  ■  0 
m 


vmach  ^ 


These  last  two  concitions  (the  variables  are  discussed  in  fef.  (1)  and 
Ref.  (2))  have  the  effect  that  there  are  no  reradiative  or  convective  losses 
at  the  surface.  The  other  quantity  needed  is  the  total  energy  required  to 
heat  the  sample  to  the  sublimation  temperature  of  the  grapMte  fibers 

12 


TIME  [s] 

Fig.  £  --  Recession  velocity  vs.  Irradiation  time  for  a  C.2F*  cm  thick 
sample  Irradiated  by  2.2  kW  cm**.  (Details  of  the  calculation  jopear  in 
the  text  above.) 


(T  «  3216*0.  This  quantity  is  calculated  in  Appendix  B;  it  is 
s  ■> 

8.2  kJ/cnT.  Upon  substituting  these  values  into  the  expressions  for 
recessional  velcity  and  burnthrough  time,  we  find 

v„  «  C.C368  cm/sec 

r ,  a  vc 

v _  *  C.C424  cm /sec 

r,max 

tg  »  6. SCC  sec. 

To  within  a  millisecond,  the  code  estimates  the  burnthrough  time  to  be 
6.9C0  sec;  the  details  of  the  burnthrough  process  can  be  seen  in  Figure  5. 


7.  COMPARISON  WITH  THE  GRIFFIS,  KASUMURA,  ANC  CKAKG  CALCULATION 
The  conditions  for  the  GHC  calculation  with  which  we  compare  are 


v„  •  C.03C6  cm/s 

r  f  avc 

V,  «  C.C361  cm/s; 

r,Bax 

these  values  are  approximate  because  we  used  the  T>3316*C  values  of  the 
surfaces  losses  to  estimate  the  energy  absorbed.  During  the  heating  process 
these  losses  would  be  somewhat  lower. 


Figure  6  can  be  compared  with  Figure  8  of  GMC.  There  is  no  substantial 
difference  between  the  two  calculations.  The  times  for  the  onset  of 
ablation  differ  by  several  tenths  of  milliseconds,  but  that  is  due  to  the 
Improved  thermal  properties  routine,  discussed  earlier.  When  I  use  the  GHC 
routine,  I  get  their  value  for  the  onset  of  ablation. 


RECESSION  VELOCITY  (cm/s) 


I 


Fig.  6  —  A  2.7S  kW/cm*  CW  Irradiation  for  comparison  with  the  equivalent 
calculation  of  Griffis,  Masumura  and  Chang.  (See  the  text  above  for  details 
of  the  Irradiation. ) 
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8.  SAMPLE  CALCULATION  FOR  A  REPETITIVELY-PULSED  IRRADIATION 

For  the  high  fluxes  that  can  occur  during  pulsed  irradiations,  a  plasma 
can  form  at  the  surface  of  the  target  material.  Cnee  it  forms,  it  absorbs 
much  of  the  incoming  laser  radiation  and  reradiates  it  over  a  broad  range  of 
wavelengths.  We  will  assume  for  this  mooel  calculation  that  half  of  the  . 
reradiated  energy  is  absorbed  by  the  graphite  epoxy.  Figure  7  is  a  log-log 
plot  of  the  front  surface  temperature  of  the  sample  as  a  function  of  the 
time  from  the  beginning  of  each  pulse;  the  results  for  the  first  two  pulses 
are  shown.  Note  the  rapidity  (less  than  C.l  us)  with  which  the  surface 
reaches  the  carbon  fiber  ablation  temperature. 


Fig.  7  —  Front  surface  temperature  during  a  repetitively-pulsed  irradiation. 
(See  the  text  above  for  details  of  the  irradiation. ) 


Figure  8  shows  the  distribution  of  the  remaining  energy  just  before  the 
start  of  the  next  pulse.  Figure  5  is  a  graph  of  the  recessional  velocity  as 
a  function  of  time  during  the  first  two  13  ^s  pulses.  We  see  that  the 
surface  recession  rate  essentially  comes  to  equilibrium  during  the  first 
Z-4  us  of  each  pulse. 

The  model  parameters  used  in  this  calculation  are 

d  *  0.1  cm 
Q  *  1.5  MW /cm2 
a  •  0 . 50 

£m  B  0,92 
vmach  *0.3 


the  next  pulse. 


ESSION  VELOCITY  [cm/s] 


TIME  [ya] 


Fig.  9  —  Recession  velocity  vs.  Irradiation  time  for  the  sample 

repetitively-pulsed  Irraaialon. 
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pulse  length  =  13  ns 
repetition  rate  =  ICC  p/s  . 

There  is  another  interesting  feature  of  these  pulsed  irradiations.  For 

our  sample  two-pulse  irradiation,  19.48  J/cm4  were  absorbed  in  the 

2 

composite  (=C.02  J/'cm  represent  the  surface  losses).  At  the  end  of  the 

-4 

two  pulses,  3.42X1C  cm  of  the  composite  had  ablated.  This  required 
54.9  kJ/cm3  X  2.42X1G"4  cm  =  18.75  J/cm2.  Thus  96  per  cent  of  the 
absorbed  energy  went  to  heat  and  ablate  the  composite,  and  only  4  per  cent 
remained  in  the  material.  This,  at  first,  surprising  result  is  due  to  the 
high  incident  fluences  and  the  low  conductivity  of  the  composite;  the 
surface  heats  up  and  ablates  much  faster  than  heat  can  flow  into  the 
interior  of  the  material. 


9.  COMPARISON  WITH  EXPERIMENTAL  REPETITIVELY-PULSED  IRRADIATIONS 

Much  of  the  data  from  experimental  repetitively-pulsed  irradiations  is 
classified,  but  we  can  use  some  of  the  data  of  Cozzens  and  Echols  [4].  The 
parameters  for  the  calculation  are 


d  a  0.2  cm 
Q  .  625  kW/cm2 
a  >  0.92 

e  3  0.92 
m 


V 

mach 
pulse  length 

repetition  rate 


0.3  . 
13  us 
ICC  p/s 


number  of  pul  ses  3  ICC 


At  this  Irradlance,  we  are  below  the  plasma  threshold,  so  we  choose  the 
absorption  coefficient  to  be  030.92;  this  is  the  same  value  used  by  GMC. 


Figures  (10M12)  show  the  temperatures  at  the  front  surface,  8th  ply 
and  15th  ply,  respectively,  for  this  16  ply  sample.  Shown  are  the 
calculated  values  and  the  experimental  values  from  Reference  4,  The  front 


surface  temper? \ures  (Fig.  (10))  were  measured  radiometrically;  the  other 
two  were  measured  with  thermocouples.  The  measured  front  surface 
temperatures,  which  were  sampled  just  before  the  onset  of  the  next  pulse, 
are  considerably  higher  than  the  Calculated  values.  They  are,  in  fact,  off 
scale  on  Fig.  (1C).  There  is,  I  think,  a  simple  explanation  for  this.  The 
code  assumes  uniform  material  at  the  surface,  which  is  slowly  eroded  by  the 
radiation.  In  practice,  there  will  be  wisps  of  graphite  fibers  sticking 
cut  that  are  not  in  good  thermal  contact  with  the  body  of  the  material. 


Fig.  10  —  Calculated  front  surface  temperature  for  a  ICO  pulse  irradiation. 
(See  the  text  above  for  details  of  the  Irradiation. ! 


TEMPERATURE  [°C) 


They  will  retain  their  heat  and  continue  glowing  much  longer  than  if  they 
were  in  good  thermal  contact  with  the  bulk  material.  Consequently,  I  think 
the  radiometer  was  picking  up  the  temperature  of  glowing  fibers,  not  the 
temperature  of  the  bulk  surface. 

The  temperatures  as  measured  at  the  8th  and  15th  plie..  are  also 
considerably  higher  than  the  calculated  values.  The  simplest  explanation 
for  this  is  that  the  conductivity  used  in  the  calculation  is  too  low.  To 
test  this  conjecture  I  arbitrarily  increased  the  low  temperature 
conductivity  by  doubling  the  values  at  10°C  and  343*C.  {The  values  1.452 
and  0.600  in  the  last  column  of  Table  A1  were  doubled.)  The  resulting 
temperature  curves  are  also  shown  on  Fig.  (11)  and  Fig.  (12).  We  see  that 
there  is  an  improvement  in  the  fit,  but  given  the  arbitrary  way  in  which  the 
conductivity  was  changed,  we  should  take  these  results  only  as  an 
Indicator. 
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TEMPERATURE  |°C] 


0  0.25  0.5  0.75  1.0 

TIME  [s] 

Fig.  12  —  Comparison  of  calculated  vs.  measured  temperature  at  the  15th  ply 
for  a  ICC  pulse  irradiation.  (See  the  text  above  for  details  of  the 
irradiation.) 

10.  . CONCLUSIONS 

The  discrepancy  between  the  calculated  and  experimental  values  is 
actually  worse  than  pictured,  for  two  reasons.  In  these  calculations,  no 
account  is  taken  of  radial  heat  flow.  The  experimental  irradiation  involves 
a  finite  sized  spot;  the  temperature  rise  behind  that  spot  will  be  less  than 
the  calculation  for  an  infinite  spot  size  would  indicate.  Additionally,  the 
sample  in  Ref.  4  had  a  protective  coating,  which  reduced  the  heat  absorbed 
into  the  material  for  several  tenths  of  seconds.  These  effects  each  lower 
the  measured  temperatures  relative  to  the  calculated  values.  In  other 
words,  all  of  the  systematic  errors  tend  to  raise  the  calculated  values  of 
interior  temperatures,  whereas  we  find  them  to  be  lower  than  the  measured 
values. 


The  most  likely  explanation  for  this  discrepancy  is  that  there  are 
errors  in  the  measured  (estimated)  [3]  values  for  the  specific  heat  and 
conductivity.  Chris  Griffis  points  out  [£]  that  the  large  peak  in  the 
specific  heat  curve  is  an  attempt  to  simulate  the  chemistry  of  the 
pyrolyzing  epoxy.  While  this  approach  may  work  reasonably  well  at  low 
fluences  and  relatively  long  time-scales,  it  may  not  on  the  short 
time-scales  associated  with  repetitively-pulsed  irradiations. 

In  any  case,  further  comparisons  of  calculation  and  experiment, 
including  any  extension  of  the  analysis  presented  here,  must  await  a  second, 
classified,  report  [6]. 
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AFP ENG IX  A:  THERMAL  PROPERTIES  CF  ThE  GRAPHITE  EPCXY  COMPOSITE 

1.  Thermal  Properties  of  Virgin  Material 

With  one  minor  nccification,  we  use  the  Menousek  anc  Monin  [2]  thermal 
properties  of  the  graphite  epoxy  composite. 


Table  A1 


Thermal  Properties  of  Graphite  Epoxy  Composite 
(Taken  from  Menousek  and  Monin  [2]) 


[NOTE: 

The  Menousek  and  Monin  version  of  the  conductivity 

runs  from  1.384  W/cm*C 

at  38*C  to  0.172 

W/cm*C  at  528*C] 

TEMPERATURE 

DENSITY 

SPECIFIC  KEAT 

CONDUCTIVITY 

TC) 

(o/ceT) 

(J/g'C) 

(W/cm*C) 

10 

1.506 

1.254 

0.C1452 

343' 

1.SC6 

2.093 

C.CG6CC 

343* 

1.EC6 

5.C24 

C.CC6C0 

510" 

1.5C6 

4.859 

C.G0172 

51C* 

1.064 

1.591 

C.GC173  , 

3316 

1.G84 

2.512 

0.CC173 

In  order  to  improve  the  numerical  stability  of  their  code,  GMC  modified 
the  Menousek  and  Nonin  properties  as  listed  in  Table  A2. 

2.  Thermal  Properties  of  Thermally  Cycled  Material 

The  large  changes  in  the  thermal  properties  between  3C0*C  and  EC0*C  are 
due  the  chemical  reactions  In  and  the  final  sublimation  of  the  epoxy.  Above 
510*C  we  assume  that  only  the  graphite  fibers  and  sot*}  residue  of  the  epoxy 
remain.  Consequently,  when  a  portion  of  the  material  that  was  heated  above 
5*0*C  cools.  Its  thermal  properties  are  different  due  to  the  loss  of  the 
epoxy.  Table  A3  lists  the  properties  for  this  case.  Any  portion  that  has 
not  reached  51G*C  Is  assumed  to  retain  all  of  Its  epoxy  and  to  still  have 
the  properties  of  the  virgin  composite.  These  properties  are  the  same 
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whether  one  starts  with  the  original  Kenousek  anc  Mcrrin  properties,  my 
modification  of  them,  or  the  £<VC  modification. 


Table  A2 

Thermal  Properties  of  Graphite  Epoxy  Composite 


(As  used  by  Griffis,  Masumura, 

and  Chang ) 

TEMPERATURE 

DENSITY 

SPECIFIC  HEAT 

CONDUCTIVITY 

CO 

(g/cm3) 

(oVg'C) 

(W/cm*C) 

10 

,1.506 

1.254 

0.01384 

329 

1.506 

2.056 

0.CC651 

357 

1.5C6 

5.0C7 

C.CC587 

496 

1.SC6 

4.902 

C.C0269 

524 

1.C84 

1.593 

C.C0228 

£66 

1.084 

1.6C7 

0.C02C6 

621 

1.064 

1.625 

C.CC187 

704 

1.084 

1.652 

C.CC179 

916 

1.C84 

‘  1.689 

0.CC173 

3316 

1.084 

2.5C8 

C.CC173 

Table  A3 

Thermal  Properties  of  Graphite  Epxoy 
That  Has  Been  Heated  above  S1C*C 


TEMPERATURE 

DENSITY 

SPECIFIC  HEAT 

CONDUCTIVITY 

Cc) 

(g/cm3) 

( J/g*C ) 

(W/cm*C) 

10 

1.C84 

1.427 

0.00173 

3316 

1.C84 

2.  ECS 

C.CC173 

26 


APPENDIX  5:  CUMULATIVE  HEAT  CAPACITY 


In  order  to  calculate  the  average  recession  velocity,  we  need  to  know 
the  energy  required  to  heat  the  graphite  epoxy  from  ambient  (2C*C)  to  the 
sublimation  temperature  of  the  graphite  (2216*0,.  Within  each  region  for 
which  the  density  and  specfic  heat  vary  linearly,  the  cumulative  heat 
capacity  is 


D(T2,T1) 


r 2 

Jji  dTp(T)Cp(T) 


(61 


(T2  "  V  UejCpj 


)  1  V  « 

1  pc 


0/.C  . 

C  pi 


+  2(52Cp2'1 


(B2 


Table  B1  gives  the  value  of  D  for  each  region  as  well  as  the  cumulative 
value  from  ambient  to  sublimation. 


Table  61 

Cumulative  Heat  Capacity  of  GraDhite  Epxoy 
(Based  on  the  Properties  in  Table  Al) 


INTERVAL 

CUMULATIVE 

TEMPERATURE 

hEAT  CAPACITY 

HEAT  CAPACITY 

■ 

D(T1-T1-1) 

0(T1#20*O 

Co 

(J/cnT) 

(J/cm3) 

20 

0 

0 

343 

821 

821 

510 

1248 

2069 

3316 

624C 

8309 

The  equivalent  value  for  the  material  after  it  has  been  heated  above  510" C 
Is  D(3316*C,2C*C).7C37  J/cm3. 


When  one  uses  the  CMC  version  of  the  thermal  properties,  one  obtains  the 
following  values  for  the  cumulative  heat  capacity. 


Table  B2 

Cumulative  Heat  Capacity  of  Graphite  Epxoy 
(Based  on  the  Thermal  Properties  in  Table  A2) 


INTERVAL 

CUMULATIVE 

TEMPERATURE 

HEAT  CAPACITY 

HEAT  CAPACITY 

D(TrTi,-1) 

0(Tr2C*C) 

Cc) 

(J/cm3) 

(J/cm3) 

20 

0 

0 

329 

770 

770 

357 

149 

919 

496 

1037 

19ES 

524 

120 

2076 

£66 

73 

2149 

621 

96 

2245 

704 

147 

2392 

816 

203 

2595 

3316 

5688 

8283 

